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INTRODUCTION TO 
PROGNOSTICS 


Motivation (1/2) 


Future aircraft systems will rely more on electrical and 
electronic components 
UAV's with all electric powertrain are increasingly being used 
for long missions 
Electrical and Electronic components have increasingly critical 
role in on-board, autonomous functions for 

Vehicle controls, communications, navigation, radar systems 


Power electronic devices such as power MOSFETs and IGBTs are 
frequently used in high-power switching circuits 


Batteries are the sole energy storage 
The integrated navigation (INAV) module combines output of the 
GPS model and inertial measurement unit. 
Assumption of new functionality increases number of faults 
with perhaps unanticipated fault modes 
We need understanding of behavior of deteriorated 
components to develop capability to anticipate failures/predict 
remaining RUL 
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“"'. Motivation (2/2) 


Images courtesy : Boeing ton, 
» — Line Replaceable Unit: ane 
Power Controller 
i A 


Lee 


Component: 
Power Transistor 


—— 


Ref: www.nasa.gov 
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Definitions 
S 


prog:nos:tic 
M-W.com — “Something that foretells” 


PHM Community —- “Estimation of the 
Remaining Useful Life of a component” 


Remaining Useful Life (RUL) - The amount 
of time a component can be expected to 
continue operating within its stated 
specifications. 

Dependent on future operating conditions 

~ Input commands 

— Environment 

_ Loads 
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Lo The Basic Idea 


Damage 

Threshold 
Damage RUL 
eT, 
EOL: 
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a Why Prognostics? 


Example: UAV Mission 
Visit waypoints to accomplish science objectives. Predict aircraft battery end of 
discharge to determine which objectives can be met. Based on prediction, plan 
optimal route. Replan if prediction changes. 

Objective #2 


Electric 
Aircraft 
\ Prognostics: 
Full discharge 
Objective #1 before mission 


completion 


Objective #3 


Objective #4 
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Why Prognostics? 


=» Prognostics can enable: 


Adopting condition-based maintenance 
strategies, instead of time-based maintenance 


Optimally scheduling maintenance 
Optimally planning for spare components 


Reconfiguring the system to avoid using the 
component before it fails 


Prolonging component life by modifying how the 
component is used (e.g., load shedding) 


Optimally plan or replan a mission 


» System operations can be optimized ina 
variety of ways 
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7 The Basic Idea Revisited 


Threshold as a Function of System State 


a a 


- Not necessarily a one-dimensional problem! = 

. ... This schematic is oversimplified! : 

System : 
State = At; = 


or 
™m 


t Time 
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Cell 
Voltage 


The Basic Idea : Batteries 
Example 


. E = End of Discharge (EOD) : 
: a 
: 
: 
: 
: 
: 
= Atrop : 
COO 2 . 
: 
: 
. 
= Voltage Threshold : 


ee 
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e 
=) 
Ls 


—_ 


The Basic Idea : Batteries 
Example 


Threshold as a 
Function of System 
State 


Future Evolution 
System of System State 
State Space 


x(te) 


What is t,? 
What is t,-t? 


What is x(t)? 
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=»—Prognostic Algorithm Categories 


Type I: Reliability Data-based 
Use population based statistical model 


These methods consider historical time to failure data which are used to 
model the failure distribution. They estimate the life of a typical 
component under nominal usage conditions. 


Ex: Weibull Analysis 


Type II: Stress-based 


Use population based fault growth model - learned from accumulated 
knowledge 


These methods also consider the environmental stresses (temperature, 
load, vibration, etc.) on the component. They estimate the life of an 
average component under specific usage conditions. 


Ex: Proportional Hazards Model 


Type III: Condition-based 
Individual component based data-driven model 


These methods also consider the measured or inferred component 
degradation. They estimate the life of a specific component under specific 
usage and degradation conditions. 


Ex: Cumulative Damage Model, Filtering and State Estimation 
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»—Data-Driven Methods 


Model is based solely on data collected from the system 


Some system knowledge may still be handy: 
What the system ‘is’ 
What the failure modes are 
What sensor information is available 
Which sensors may contain indicators of fault progression (and 
how those signals may ‘grow’) 
General steps: 
Gather what information you can (if any) 
Determine which sensors give good trends 
Process the data to “clean it up” - try to get nice, monotonic 
trends 
Determine threshold(s) either from experience (data) or 
requirements 
Use the model to predict RUL 
Regression / trending 
Mapping (e.g., using a neural network) 
Statistics 
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»—Data-Driven Methods 


Pros 
Easy and Fast to implement 
Several off-the-shelf packages are available for data mining 
May identify relationships that were not previously considered 
Can consider all relationships without prejudice 


Cons 


Requires lots of data and a “balanced” approach 
Most of the time, lots of run-to-failure data are not available 
High risk of “over-learning” the data 
Conversely, there’s also a risk of “over-generalizing” 
Results may be counter- (or even un-)intuitive 
Correlation does not always imply causality! 
Can be computationally intensive, both for analysis and implementation 


Example techniques 
Regression analysis 
Neural Networks (NN) 
Bayesian updates 
Relevance vector machines (RVM) 
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==—Physics-Based Methods 


Description of a system’s underlying physics using 
Suitable representation 


Some examples: 


Model derived from “First Principles” 
» Encapsulate fundamental laws of physics 
PDEs 
Euler-Lagrange Equations 
Empirical model chosen based on an understanding of the 
dynamics of a system 
Lumped Parameter Model 
Classical 1st (or higher) order response curves 
Mappings of stressors onto damage accumulation 
Finite Element Model 
High-fidelity Simulation Model 


Something in the model correlates to the failure 
mode(s) of interest 
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»——Physics-Based Models 


Pros 
Results tend to be intuitive 
Based on modeled phenomenon 
And when they’re not, they’re still instructive (e.g., identifying needs for more 
fidelity or unmodeled effects) 


Models can be reused 
Tuning of parameters can be used to account for differences in design 


If incorporated early enough in the design process, can drive sensor 
requirements (adding or removing) 
Computationally efficient to implement 


Cons 
Model development requires a thorough understanding of the system 


High-fidelity models can be computationally intensive 


Examples 
Paris-Erdogan Crack Growth Model 
Taylor tool wear model 
Corrosion model 
Abrasion model 
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INTRODUCTION TO 
MODEL-BASED 
PROGNOSTICS 


=. Model-based prognostics (1/2) 


RUL Pa Failure 
se Computation Threshold 
Do {a, B} 
RUL (ty) 
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State vector includes 
dynamics of the 
degradation process 


It might include nominal 
operation dynamics 


EOL defined at time in 
which performance 
variable cross failure 
threshold 

Failure threshold could 
be crisp or also a random 
variable 


", Model-based prognostics (2/2) 


=» Tracking of health ~~ ees 


-- 3€-- Measured 
state based on || Prete ae [t,=116] 
measurements sf - | : , 


110 120 130 140 150 160 170 180 190 


=» Forecasting of 20, 


health state until ai| 
failure threshold is = ,.-# 
crossed 20 
= Compute RULas || | 
function of EOL eke 
defined at time 20 
failure threshold is * . 
crossed oe a 


5 i i i J 
110 120 130 140 150 160 170 180 190 
Aging time (hr) 
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“. Methodology 


Accelerated 
Aging 


Training 
Trajectories 


Lp = AXp—-1 + Bup_1 + We-1 


Ye = Hr, + VE 


Degradation 
Modeling 


Parameter 
Estimation 


Health State RUL 
Estimation Estimation 


Test 
Trajectory 


State-space 
Representation 


Dynamic 
System 
Realization 


Kalman Health State {@(tp+1)s +++ €(tptn)f 
Filter Forecasting 


RUL ns Failure 
Computation Threshold 


RUL(ty) 


21 
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RESEARCH APPROACH 
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ial 


High level research efforts 


Prognostics models and algorithms 
» Identification of precursors of failure for MOSFETs under different failure 
mechanism conditions 
Identification of precursors of failure for different IGBT technologies 
(CALCE) 
Modeling of degradation process MOSFETs 
» Development of prognostics algorithms 
Prognostics for output capacitor in power supplies (Vanderbilt 
University) 
» Electrical overstress and thermal overstress 
» Development of prognostics algorithms 
Accelerated Life Testing 
» Thermal overstress aging of MOSFETs and IGBTs 
» Electrical overstress aging testbed MOSFETs 
» Electrical overstress aging testbed for Capacitors 
Effects of lightning events of MOSFETS (LaRC) 
Battery Degradation and ageing ( ARC — LaRC) 


Ageing Effecting on ESC’s ( ARC — LaRC) 
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Research Approach 


Kolin luha(er-haceyame)mr-|i|0)q=m anvele(=t-mr-|avemual=)| quasi t-la(o)atsjal|emuemual-lienel-laulellt-le 
g-J1iUlasmnalsvelat-lalisiaals 


Identification of precursors of failure which play an essential role in the prediction of 
ig=jaat=lialiale miK= 


DY=\V{=) (0) 9) a nk=) a] ee) mr=1XeC>) (=) e-1u=(omr- (6) ale mu =i-)n ol -1el-mu al-lemr-lell /ju-]K-m lal omy 40) (ele-lule)ame)mellir-la-lale 
failure mechanisms and aid the understanding of damage progression 


Development of degradation models based on the physics of the device and the failure 
inal=eiat-lalisiaals 


DY=\V{=1 Ko) o)aal=1alemedmma=iaat=lialiale Mm lixowe) avoid (eamr-] (ele) siaalan\-mual-leme-] <i ipinem-[eoelelalmmnal= 
different sources of uncertainty while leveraging physics-based degradation models 
iu at=] Mere) al-j[e(~)a>fnn0] G0] q-me)98\~)a-hule)al-)mr-] elem =lanvdiae)alan(=lale-]meae)alelinte) als 
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Prognostics Algorithm Maturation 
through Validation Experiments 


DY =¥-4e-(eF-1alelam \VAlelel=)| 
Development 


Prognostics Algorithm 
Update 


Validation experiments 


ime) ganlel -hulelameymanteye(=)m-lare| 

Fl ieXolaidalantcminalo)genyclealsales WE) ifekclulelaMmelarel=lane) ekt-1avccre| 

Fale mel=t:4e-\er-|nle lam oceans experiments 
hypothesis 
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Prognostics Algorithm Maturation 
through Validation Experiments 


Aging under therMd! een 


and electric = f 
Aging under thermal aii, av 


TRL4 
Aging under Oa stresses 


TRL 3 - 


i ¢) — ( 2ene Veo— Ve (#) 
a o ae (a) Gee, 
= av 


ESR(t) = 5 (ep do Pz) (ge re | 


First Principles 


Oo e ESR(t) = 1 (pe dc Pr) ( Jeo t We ) Model for multiple failure 
Ss e 


Veo—Ve (t) modes 
First Principles 
; Model 
= )= 4 +B 
Physics inspired 
empirical 
Model 
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ELECTRIC VEHICLE 
POWERTRAIN 
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System Level Prognostics 


=» Component Level Prognostics 
—- System Level Prognostics 
» Batteries 


» Power Conditioning Circuit — 
Capacitors, MOSFETs 


Electronic Speed Controllers ( 
ESC) - MOSFETs 


» BLDC 
=» Study Cascading faults 


=» Effects of component level 
aging/degradation on system 
performance 
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7 Hardware in Test Loop 
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ACCELERATED AGING AS 
A PROGNOSTICS 
RESEARCH TOOL 


Accelerated Aging 


Traditionally used to assess the reliability of products with 
expected lifetimes in the order of thousands of hours 


in a considerably shorter amount of time 
Provides opportunities for the development and validation of 
prognostic algorithms 
Such experiments are invaluable since run-to-failure data for 
prognostics is rarely or never available 
Unlike reliability studies, prognostics is concerned not only 
with time to failure of devices but with the degradation 
process leading to an irreversible failure 
This requires in-situ measurements of key output variables and 
observable parameters in the accelerated aging process with the 
associated time information 
Thermal, electrical and mechanical overstresses are commonly 
used for accelerated aging tests of electronics 
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Example: Electrical 
overstress aging of Power 
Transistors 
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Accelerate aging strategy (1/3) 


=» The main strategy is the 
application of electrical overstress 


fixed junction temperature in order to avoid 
thermal cycles 


avoid package related failures 


=» Accelerated test conditions are achieved 
by electrical operation regime of the 
devices at temperatures within the 
range below maximum ratings and 
above the room temperatures. 
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a Accelerate aging strategy (2/2) 


=» The highest acceleration Simulated I-V characteristics and 
factor for aging can be instability boundary at 300° K for 
achieved in the proximity of power MOSFET. 
the SOA boundary 


=» Instability points represent 
the critical voltages and 
currents limiting the SOA 


=» An electrical regime close to 
the SOA boundary serves as 
the accelerator factor 
(stressor) and it is expected 
to reduce the life of the 
device 

=» The safe operation area 


boundary shifts closer to 


Rosi 0 20 40 60 80 §=10 
the origin as the DRAIN-SOURCE VOLTAGE (V) 
temperature increases 
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IN 


1) 


Bie Ee ieee ek ROSE BA GOS 5, OS ak Re eee Ee a eR Gk RBA NN 2 eee ae eS & 


— 


oO 


DRAIN CURRENT x100000 (A/um) 
NO 


34 


Aging system description (1/3) 


=» Three main components in terms of hardware 
Electrical operation unit of the device 


» custom made printed circuit boards for the 
instrumentation circuitry and gate drivers 
» commercially available power supplies and function 
generator to control the operation of the DUT 
An in-situ measurement unit of key electrical and 
thermal parameters 
» commercially available measurement and data 
acquisition for slow and high speed measurements 
Thermal block section for monitoring and control of the 
temperature 
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a 


Aging system description (2/3) 


Thermal block for measurement and control of 
device temperature 


MO] 29H 


Thermo-electric Unit Thermocouple 
Modules 


Heat Sink with Fan 
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Aging system description (3/3) 


Experiment on power MOSFET (1/2) 


IRF520Npbf power MOSFET 
» TO220 package,100V/9A. 
Electrical overstress used as acceleration 
factor. High potential at the gate 
Vgs=50V, Vgs rating is 20V max. 
Vds=2.4V with a 0.2 ohm load. 
Temperatures kept below maximum rating 
T;max=175°C 
Objective is to induce failure mechanism on 
the gate structure 
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Experiment on power MOSFET (2/2) 


Degradation process as observed on 
threshold voltage (V,,) 


Threshold Voltage Measurements 


2.5 |. 


nD 
] 


Current I (A) 


0.5 |. 


—t— Pristine 
—cs— Aging1 
—+— Aging2 
©?) —e— Aging3 || 
| Y | —2— Aging4 
4 Aging5 
y~ | |——Agings 
7— Aging7 


“— Aging8 
—+— Aging9 
wena ageree: Aging10 || 


Vii, shifts right as a 
result of degradation 


Fa * Devices lost gate 
appt control and failure is 

ri ive oe ar “iFFeve rsi b l e 

Voltage Vos=“as (V) 
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Example: Electrical 
overstress aging of 
Electrolytic Capacitors 
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Accelerated aging system 


= Allows for the understanding of the effects of 
failure mechanisms, and the identification of 
leading indicators of failure essential for the 
development of physics-based degradation 
models and RUL prediction 

= Electrolytic capacitor 2200uF, 10V and 1A 

= Electrical overstress >200 hr 


Square signal at 200 mHz with 12V amplitude and 
100 ohm load 
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Electrical Overstress Aging 


* System 


Square wave amplified 
signal 
Power Supply | | R. 
= 


Signal Amplifier 
hardware ‘ 
Agilent Signal | 


Generator Vo 


Input Square wave 


C 
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"" 


-Im(Z)/Ohm 


Degradation observed on EIS 
measurements 


0.06 
-Im({Z) vs. Re{Z) 
0.05 oe Pristene.mpt 
, =< 71hours.mpt 
se 161Hours.mpt 
0.04 gf <= 194Hours.mpt # 
0.03 
0.02 
0.01 
0 
-0.01 
-0.02 
-0.03 
-0.04 


Re(Z)/Ohm 


CASE STUDY I: 
PROGNOSTICS OF 
ELECTROLYTIC 
CAPACITORS 


a Case Study: Avionics System 


=» Integrated Avionics systems consists of: 
Global Positioning System (GPS) module 


Integrated navigation (INAV) module combines 
output of the GPS model and Inertial measurement 
unit 

» Power Supply module 


Fault Insertion Points 


GPS: Global Positioning System 
INAV: Integrated Navigation 
PSM: Power Supply Module 
IMU: Inertial Measuremen t Unit 
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a Methodology 


Accelerated 
Aging 


Training 
Trajectories 


Degradation 
Modeling 


LE = ATR-1 + Bug—i + We-1 


Test Uk = va kigs —+- Uk 


Trajectory 


State-space 
Representation 


Dynamic 
System 
Realization 


{E(tp41); as ,&(tp+n)} 


RUL <- Failure 
Computation Threshold 


Health State RUL 
Estimation Estimation 


RU L(tp) 
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Accelerated Aging and 
Precursors of Failure 
Features 
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Electrical Overstress Aging 
System 


Square wave amplified 
signal 
Power Supply | | R 
| L 
Signal Amplifier 
hardware ‘ 


Agilent Signal C 
Generator Vo Vi 


Input Square wave 
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-Im(Z)/Ohm 


Degradation observed on EIS 
measurements 


¢ -Im{Z) vs. Re{Z) 
ane j = Pristene.mpt 
, —= 71hours.mpt 
| =e 161Hours. mpt 
0.04 f <= 194Hours.mpt # 
0.03 Pi 
0.02 $ 
0.01 ff 
0 fF 
0.01 } 
-0.02 q 
-0.03 
we 
-0.04 
0.05 O41 015 0.2 
Re(Z)/Ohm 
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Degradation on lumped 
parameter model 


ESR C and ESR are estimated 
— /-—W—« from EIS measurements 


C 
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ote] | 
Capacitor 


Methodology 


Reduced 
order 


Degradati 
fo) am aaleyersy| 
(time) 


Health 


state 
estimation 


RUL 
prediction 
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Empirical degradation model 


Based on observed degradation from 
Capacitance parameter 

Using training capacitor data to 
estimate degradation model parameters 
Assumed exponential model based on 
Capacitance loss 

Parameter estimation with least- 
Squared regression 


C.=e “+f ! fA 
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a. Degradation model results 


0.0162 -0.8398 
cee MF | (0.0160, 0.0164) | (-1.1373, -0.5423) | 1°77® 
0.0162 -0.8287 


0.0161 -0.8217 

m-8.45.46 | ois, 162 | ca.zs 050m | 18860 
0.0162 -0.7847 

mI#*, #0 | (0.0161, 0.0164) | 1.1134, -0.4560) | 714! 
0.0169 “1.0049 

= #I® | 0.0167, 0.0170) | (1.2646, -0.7453) | 77°! 


The optimal parameter presented along the 95% confidence 
interval. 


e The residuals are modeled as a normally distributed random 
variable with zero mean and variance 
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Prognostics algorithm 


Implementation of prognostics algorithm with 
Kalman filter 


Capacitance loss considered as state variable 


EIS measurements and lumped parameter model 
used to obtained measured capacitance loss values 


Empirical degradation model used to generate the 
state transition equation 


Use one Capacitor for testing and the rest for 
model parameter estimation (leave on out test) 


Failure threshold of 20% drop on capacitance 
based on MIL-C-62F 
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a Kalman filter implementation 


=» State transition equation 


» State-space model for 
derived from degradation filter implementation 
model at, 


Ci=e “+p 
Ul 
0C 


—=aC -ap 


Cr. = ArCx -1+ Biut+v 


at ==>) ye =hCk+w,where 
Ci—Cr- ar 


Ax = (1+ As), 
B. =-ap,, 


=AC:-a—-C 
At P 


Cis ( + OA:)C: - At — apa: 
Ce = (1+ GAk)Ce - 1 — Ak 


h=lu=1. 
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Prediction mode 


Assumed measurements are not 
available at some point in time 


Filter used in forecasting mode to 
predict future states 


Predictions done at 1 hr. intervals 


State transition equation used to 
propagate state (n: number of 
prediction steps, /: last measurement at 
[ 7 n-\ . 

; C4, =A°C, + SA'B 


l+n 
i=0 
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| Tracking and forecasting (Cap. 
7 #6) 


-«-- Measured 
15;| © Filtered 
Predicted 


— eee’ (alpha = 0.3, lambda = 0.5) 


110 120 130 140 150 160 170 180 190 


110 120 130 140 150 160 170 180 190 


110 120 130 140 150 160 170 180 190 


10} a ee 


5 La ! 1 ! nl 1 
110 120 130 140 150 160 170 180 190 
Aging time (hr) 


57 
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7 Relative Accuracy 


RUL* — RUL 
RU L* 


94 | 85.6 | 90 | 789 | 94.8 | 94.2 | 90_| 


Tie | 6 [98.1 [765 [98 | 96.2 | 96 


“16.6 [87.5 | 808 | 87.34 


Fiase [seawater ao 
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RA = 100{1- 


58 


End of Case Study I: 


QUESTIONS? 
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CASE STUDY II: 
PROGNOSTICS OF POWER 
TRANSISTORS 


Modeling for Power MOSFET under 


- electrical overstress 


om GatSource 


- Two-transistor model is shown to be 
a good candidate for a degradation 
model for model-based prognostics. 


- The model parameters K, and W1 
could be varied as the device 
degrades as a function of usage time, 
loading and environmental 
conditions. 


- Parameter W1 defines the area of the 
healthy transistors, the lower this 
area, the larger the degradation in Vop  Wy=0.9 W,=0.1 
the two-transistor model. In addition, 
parameter K serves as a scaling 
factor for the thermal resistance of 
the degraded transistors, the larger 
this factor, the larger the 
degradation in the model. 
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Precursor of Failure 


As case temperature increases, 
ON-resistance increases ' 


Recrcpy ¥S- Temperature: Device 9 


This relationship shifts as the eli: 
degradation of the device 


i@eeegees 


increases Ee 

For a degraded state, ON- * sl 

resistance will be higher at any). 

given case temperature aap" 

This is consistent with the die- °*" ms 


attach damage since it results 
On increased junction 
temperature operation 

This plot can be used directly 
for fault detection and 
diagnostics of the die-attach 
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a Degradation process data 


Normalized ON-state resistance (ARpsion)) 
and filtered trajectory for device #36 


Normalized ON-state resistance (ARpgon)) 


—— ROS on and filtered trajectory for device #36 


Filtered 


Cases #08, #09, #11, #12 
and #14 are used for algorithm 
development purposes. 

Case #36 is used to test the 
algorithms. 
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a Empirical Degradation Model 


=» An empirical degradation model was selected for the 
model-based algorithms 


=» Exponential based function to capture degradation 
process 


=» Two parameters in the model which will be estimated 


==> ARpson) = @(e =I) 
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Prediction of Remaining Life 
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RUL Prediction Methodology 
Considerations 


A single feature is used to assess the health state of the 
device (ARpscon)) 

It is assumed that the die-attached failure mechanism is the 
only active degradation during the accelerated aging 
experiment 


Furthermore, ARpsony accounts for the degradation 
progression from nominal condition through failure 


Periodic measurements with fixed sampling rate are available 
for ARpsvon) 


A crisp failure threshold of 0.05 increase in ARpsony is used 


The prognostics algorithm will make a prediction of the 
remaining useful life at time ¢,, using all the measurements up 
to this point either to estimate the health state at time ¢, ina 
regression framework or in a Bayesian state tracking 
framework 


It is also assumed that the future load conditions do not vary 
Significantly from past load conditions 


VII International School on Fault Diagnosis of Complex Systems (Tarrasa, July 2017) 66 


a RUL Prediction Algorithms 


=» Gaussian Process Regression 


Algorithm development cases used to select covariance matrix 
structure and values 


» Extended Kalman filter 
Empirical degradation model 


State variable: Normalized ON-resistance and degradation model 
parameters 


Arbitrary values for measurement and process noise variance 
= Particle filter 
Empirical degradation model 


State variable: Normalized ON-resistance, degradation model 
parameters 


Exponential growth model used for degradation model parameters 
Arbitrary values for measurement and process noise variance 
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RUL 


120 ., 


100 |. 


Ok 


RUL estimation results 
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End of Case Study II: 


QUESTIONS? 
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CASE STUDY III: 
PHYSICS-BASED 
PROGNOSTICS OF 
CAPACITORS 


a. Capacitor Structure 


=» An aluminum electrolytic capacitor, consists of 
» Cathode aluminum foil, 
Electrolytic paper, electrolyte 
» Aluminum oxide layer on the anode foil surface, which acts as 
the dielectric. 
Equivalent series resistance (ESR) and capacitance(C) are 


A 


electrical parameters that define capacitor health sco. eidieinae 


foil 


Connecting Lead : ; 
g dielectric Layer 


amet 7) PSE INN G 
: Ye Uj 
highly etched La a i 

4 i, 


aluminum foil 


Anode Foil 


ay Foil 


Separator Paper 


leakage current 


Al2O3 — 
electrochemical 
oxide 
layer(forming) 


electrolyte paper Al,O3 — oxide 
(spacer) layer(natural) 


Physical Structure Internal Structure 
Ref :http://en.wikipedia.org/wiki/File: ElectrolyticCapacitorDisassembled.jpg 


Open Structure 
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a Degradation Mechanisms 


Electrolyte Evaporation 


, Aging i tle ; Over Voltage Stress 
dielectric material 


Excess Ripple Current 


Degradation in 
Anode foil 
Charging\Discharging Cycles 
alk iki Over Voltage Stress 
High Ambient Temperature : 
Excess Ripple Current 


Prolonged Use _ 
Increase in internal Temperature Chaiging\Discharging Cycles 


Electrical Stress 


Prolonged Use -Nominal Degradation i iL 
Over Voltage Stress 


Decrease in capacitance 
F Over Voltage Stress Increase in ESR 
Excess Ripple Current 


Charging\Discharging Cycles Excess Ripple Current 
High Ambient Temperature 


ous Charging\Discharging Cycles 
Increase in internal Temperature 


Degradation in 
Prolonged Use -Nominal Degradation Cathode {oll 
Degradation of Oxide Film 


Degradation Causes\ Mechanisms Failure Modes 


72 


VII International School on Fault Diagnosis of Complex Systems (Tarrasa, July 2017) 


a Capacitor Degradation Model 


Pristine Capacitor a 
= dielectric Layer foil 
| ES R F highly etched leakage current 
— | e | /—-w—s Deg ra d. aluminum foil 
C C M 1 sletootemia | 
layer(forming) 
Id €a N Oo n = ” ; electrolyte paper Al,03 — oxide 
ei I d e a | (spacer) layer(natural) 
lectrolyte volume V, maximum Th ' 
Capacitance Value maximum erma 


tress 


Electrolyte degradation + 
Decrease in (A,) + 
crystallization +oxide layer 
breakdow 


Coxide_layer 


Avg. surface area decreases 
(A,) + oxide layer 


2maQ 


0.1 


lagnosis of Complex Systems (Tarrasa, J . Fe) (ohm) 


Sl 
i 


Empirical Model with static 
parameters 


This empirical model represents an approximation of 
lumped parameter model 
E, : C(t) =e“ + B, 
e aand 8 are degradation model parameters estimated from the 
experimental data. 
The following system structure is used in the 
implementation of the filtering and the prediction using the 


(1 + Ag), 


Kalman filter. B, = ofA, 


h I 
t=1 


The state variable (x,) at aging time (t, ) is the percentage 


loss in Capacitance. 
Process noise was estimated from the model regression for the empirical model 
Measurement noise was estimated from the EIS measurements 
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Degradation Model: Electrical 
Circuit Equivalent 


Cc 


". 


Al2O3 — oxide 
layer(natural) 


electrolyte paper 
(spacer) 


a 1mQ 1mQ >» 2mQ 
O ; 
= 


Roao2 10K Ce Roo210K “4 Ro 


Coxide_layer \ 
» 


Granularity of 
Degradation 
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Capacitance Degradation Model 
"" 


=» Decrease in electrolyte volume : 
(1) 


=» Capacitance (C) ): Physics-Based Model: 


a (2) 


=» Electrolyte evaporation dominant degradation phenomenon 


» First principles: Capacitance degradation as a function of 
electrolyte loss 


Z2€REO Veo — Ve(E) 
D, : C(t) = | —— } (——_—-}, 3 
\ ( de ( Jeo t We ) ( ) 
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Capacitance Degradation 
Model 


Oxide breakdown observed - experimental data 


The breakdown factor is exp. function of electrolyte 
evaporation 


Coxct) = exp f(V., - Vecty) 


Updated in capacitance degradation model : 
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7 Dynamic Model of Capacitance 


(4) 


(9) 
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7 Dynamic Model of Capacitance 


(6) 


(7) 
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a Dynamic Model of ESR 


=» Decrease in electrolyte volume : 


Vel(t) = Vix tt Agden 


» ESR 


» Based on mechanical structure and electrochemistry. 
» With changes in R, (electrolyte resistance ) 


Ps: shen ~ Esme ~ (sr Ped caw) (8) 
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=3a—_ ON-LINE——_> t+? OF F-LINE———> 


Process Flow 


Monitor degradation 
behavior 


Physics-based / 
Empirical 
Degradation models 


Offline Parameter 


Estimation 


Health State 


Tracking Health 
State /Parameter 
Estimation 


Compute 
RUL 
Performance 
Assesment 


RA = 100 


D4: Crzi = Cr ( 


2¢ RE ow, Ask oCbhk 


(C F ; Jeot i)> Cbk(i) ) 


(C) 


J) 
de. 


RU L* — RUL 


RU L* 
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Unscented Kalman Filter for 
" State Estimation 


2€REQWeAsJeoCbk 
D4 ; Cr+1 = Ck _~ (Set At 
Cc 


=» Derived physics-based degradation model 


=» The following system structure is implemented for state 
estimation i 


=» The state variable (x,) is the current health state at 
aging time (t,) 
Process noise was estimated from the model regression for the empirical model 


Measurement noise was estimated from the EIS measurements 
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Electrolyte Volume Estimation 
for TOS Experiment 


Parameter | xX | x | $.D | Cl 

6,(mm*) 523.6112 | 523.6113 0.0026 (523.6098, 523.6127] 

62(mm?* /t) 0.0161 0.0161 ] 1.8748 x 10°" | [0.01614, 0.01611} 

63(mm/t?) 3.8077 x 10" | 3.8072 x 10” | 6.9373 x 10” [0.3769 x 107°, 0.3846 x 107°) 
RMSE 26.2232 [ 26.2277 | 0.0483 (26.1965, 26.2500) 

RMSPE 0.8589 | 0.8591 |. 0.0016 10.8580, 0.8598 


Summary for Linear Regression Electrolyte Degradation Model 


T 


« Volume Decrease ( Cap. #1 -15) 
—e-— Estimated volume (fit) 


Volume (mn? ) 
Residuals 


Residuals | 


500 = 1000 1500 2000 2500 3000 3500 500 1000 1500 2000 2500 3000 
Aging Time ( Hours ) 2 
Time ( Hours ) 


Volume estimation for all capacitors Volume estimation Error 
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RA Results — Discussion 
“/ EOS Experiment 


Capacitance - Over RA summary for model ¢, 


Aging | RA, 
Time 


Ej ° C(t) = 


100 
Time 


Capacitance - Over RA summary for model D4 


a=0.3, B=0.5 


2 We A. leo “OK 
Da: Crti = Cr — (Fee | At 


d2, 
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Time (s) 


Tracking 


Capacitance q F) 


Alpha Lambda 


Capacitance (,, F) 


RUL and Validation — EOS - 
Experiment — Capacitance 
Degradation Model 4 


| —— Measured Predicted | 


—s— measured data 
+ filter data 


fa» ae 


ae 


2 
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oe 
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r L 
180 =. 200 
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Predictions at different aging time 
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RUL and Validation -— EOS -Experiment 
" - ESR Degradation Model ?: 


Cap #2 


—s— measured data 
+ filter data 


a i 
~ * 


L L L L L L t 
80 100 120 140 160 180 200 
Aging Time ( Hours) 
Output Error - Cap #2 


Tracking 


r 


[ ie 
140 160 180 


t t t t t t t t t 
40 60 80 100 120 140 160 180 200 
Aging Time ( Hours) 


a=0.3, B=0.5 t r t r r 
100 120 140 160 180 200 


Predictions at different aging time 


t r r t 
140 160 180 200 
I =171 


RUL (hours) 


Alpha Lambda 


[ c c [ c c Cc 
80 100 120 140 160 180 200 
Aging Time ( Hours) 


80 100 120 
Aging Time (Hours) 
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Alpha Lambda 


RUL and Validation — TOS - 
"s Experiment - Capacitance 


3 
x 10 
pies 


3 & Measured Predicted | 
x 10 


—=— measured data | 
+ filter data 


Pre an 
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8h SE bbb be teeice iy ie 


ieee t, = 87 
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End of Case Study III: 


QUESTIONS? 
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CASE STUDY IV: 
PROGNOSTICS OF LI-ION 
BATTERIES 


a Battery Modeling 


— Equivalent Circuit Empirical Models 
= Most common approach 


» Various model complexities used 
. Difficulty) in incorporating agiNc (p) C 


R_ ohmic |_batt R_ohmic 


OCV(SoC) | U OCV(SoC) 


(c) Es el Gs C_| 


R_ohmic 
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«Battery Model — Tuned using Lab Data 


NW WW 


=» An equivalent circuit battery model is ‘- 
used to represent the battery terminal . 
voltage as a function of current and the 
charge stored in 3 capacitive elements = 
X = [4b ep gcs|" — V Measured 
my Fitted 
0 0 0 4 
r= | 0 = 0 r+] 1 Jité 
0 OO ~-RE 
y=v=[2 _— | -2|-* 16 
Cp Cep Gs 0 10 20 30 40 50 
, , Time (h 
= Two laboratory loading experiments are a asia 
used to fit the following parameterization SS — V Measured 
Di “~~ B3 V Modeled === 'V Fitted 


coefficients 
SOC =1-— max — Yb 


Y 
C mar 


( 'b =¢ i vo + € Y, ‘bl “SOC +¢ Y, *b2 ‘SOC* +¢ 3° SOC : 


Cop = Copo + Cepi + XP (Cope (1 — SOC)) 


Rep = Repo + Repi : exp (Repa (1 — SOC)) 


2 3 
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a Battery Modeling 


— Electrochemical Models vs. Empirical Models 
» Battery physics models enable more direct representation of age-related changes in 
battery dynamics than empirical models 
= Typically have a higher computational cost and more unknown parameters 


+— 


i(t) Load 


Charge 
Reduction at pos. electrode: Oxidation at pos. electrode: 
a ad + ne > LiCoO, > Li,_,CoO, + nLit + 
LiCoO, ne- 
key al - — saree Reduction at neg. electrode: 
I, nly + ne + nLit + ne +C>Li,C 
Current flows + to - Current flows - to + ° 
eect OTP intermayjonsi choo on Fault Diagnosis of Complex systenta CSO SLAPS 1) to - 


nrnc flraivs Ltn — 


a Electrochemical Li-ion Model 


=» Lumped-parameter, ordinary differential equations 
=» Capture voltage contributions from different sources 
» Equilibrium potential >Nernst equation with Redlich-Kister 


expansion 
Concentration overpotential > split electrodes into surface and 
bulk control volumes aie: 


Surface overpotential > 
Butler-Volmer equation 
applied at surface layers 


Ohmic overpotential > “ 
Constant lumped resistance —_ "smn cA 
accounting for current ——= —— 
collector resistances, 
electrolyte resistance, 
solid-phase ohmic resistances 
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Battery Aging 


Simulated 


=» Contributions from both decrease in —~., 
mobile Li ions (lost due to side reactions 
related to aging) and increase in internal 
resistance 
Modeled with decrease in “g”*” parameter, 
used to compute mole 1 a 


Modeled with increase in “Rk,” parameter 
capturing lumped resistances 


Cycle 16 
Cycle 26 
Cycle 36 
Cycle 46 
Cycle 56 
Cycle 66 
Cycle 76 
Cycle 86 
Cycle 96 
Cycle 106 
Cycle 116 
Cycle 126 
— Cycle 136 
——— Cycle 146 
———— Cycle 156 
Cycle 166 


254 : IN : : Cycle 176 
2000 4000 6000 8000 Cycle 186 


Current (A) 
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i Measured 


Voltage (V) 


"s 


Subscale electric 
aircraft operated at 
NASA Langley 
Research Center 


Powered by four sets 
of Li-polymer 
batteries 

Estimate SOC online 
and provide EOD and 
remaining flight time 
predictions for 
ground-based pilots 


Edge 540-T 


Powertrain 


5 Ea + 
Ea 


582P + . 4 
LiPo 
t12 


| M2 Propeller 
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a. Edge UAV Use Case 


=» Piloted and autonomous 
missions, visiting waypoints 

=» Require 2-minute warning for 
EOD so pilot/autopilot has 
sufficient time to land safely 


om _ 
This answer depends on battery age Arcee »* — _ 

» Need to track both current level of x mae 
charge and current battery age Objective #1 ‘y 
Based on current battery state, current _ 6 
battery age, and expected future usage, _ oo eeu - 
can predict EOD and correctly issue 2- Obiective oe X Eu” 
minute warning <-=---" 

Runway 


VII International School on Fault Diagnosis of Complex Systems (Tarrasa, July 2017) 


- 


Predicted Time Remaining (min) 


¢ The min, max and median 


20; 


nn 
T 


Measured and predicted 
battery current, voltage 
and SOC different time 
steps 


predictions are plotted 
from each sample time 
until the predicated SOC 
reaches 30% 


Median Predicted 

aseese Time Remaining (Ground Truth) 
—— Range Predicted 
fey. 
I: 
a * ts, 
he 

0 5 10 45 


Time (min) 


Predication over Flight Plan 


—— Median Predicted ='='= Min/Max Predicted ——— Measured| | Rp On 


Prediction at t = 0.02 min Prediction at t = 2.62 min Prediction at t= 5.15 min 


age (V) Pack Current (A) 


Lowest SOC (%) Pack Volt: 


0 5 10 15 0 5 10 15 0 5 10 15 


Time (min) Time (min) Time (min) 


Predictions for remaining flight time for 
entire flight plan 

Overestimate till parasitic load is injected 
Once the parasitic load is detected the 
remaining flying time time prediction shifts 
down. 


Ref : E. Hogge et al, “Verification of a Remaining Flying Time Prediction System for Small Electric Aircraft”, 
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= 
= 


uv 


Performance Requirements 


=» Accuracy requirements for the two minute warning were specified as: 


» The prognostic algorithm shall raise an alarm no later than two minutes 
before the lowest battery SOC estimate falls below 30% for at least 90% of 
verification trial runs. 

» The prognostic algorithm shall raise an alarm no earlier than three minutes 
before the lowest battery SOC estimate falls below 30% for at least 90% of 
verification trial runs. 

» Verification trial statistics must be computed using at least 20 
experimental runs 


8 T T T 25 i * i * 
c Median + u 25%-75% 9%-91% - : & : 
6F i I * * I 
a i : e ! 
er — La 3.4 
S 4, —_ T ii f=] i] * * i] 
is ee ee -: ae sy ri a : @! [=== Thresholds 
5 2+ 1 “i 1 ot bi : 1 _ 7 & 15 rl * * 1 * = = With parasitic load 
¥g Ae T 5 ; * : * = = Without parasitic load 
—E OF L 1 a u * os * = Out of Bounds 
= a 1 5 1 * 1 
iz +L af Port > 10 1 & 1 
= 2 i i+ L i r * | 
- I 
9 a - - I = 1 i” ‘ mean = 2.4e+00 
i] | 47) ¢ . 
al 2% i a 5 ? p (ee ase 
i] i] 
i] i] 
-6 — L . Ll I * I 1 1 1 
20 25 30 35 40 4 3.5 3 2.5 2 1.5 1 0.5 0 
BatteryID Remaining Useful Life (min) 
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End of Case Study IV: 


QUESTIONS? 
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Data Sets Available for 
Download 


eee |) | /Tti. arc.nasa. gS as ‘prognostic-data- 
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ee eee ee ee , 
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CLOSING REMARKS 
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Remarks (1/2) 


=» Electrical and Electronics PHM Maturity - 
scientific and engineering challenges 
=» Research approach challenges 
How to balance lack of knowledge of the 
system vs own expertise on particular PHM 
tools 
- Data-driven or model-based? 

» Data is always needed but more 
important, information about 
degradation/aging processes is key 

» Experiments and field data 
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Remarks (2/2) 


Aging systems as a research tool 


Value in terms of exploration of precursors of failure 
and their measurements is evident 


Still an open question on how degradation models 
and algorithms are translated to the real usage 
timescale 

In the use of physics 
It should be embraced 


Validate models and algorithms with data from 
lab experiments and fielded systems 


A success in developing PHM methodologies in 
an real usage application will require the right 
team 
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